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I. INTRODUCTION 

We define the optical Poynting singularity as an iso- 
lated point, line or surface at which the Poynting vector 
of electromagnetic beam vanishes [l(. This singularity is 
understood in the sense of indeterminate direction of the 
Poynting vector @. In Refs. 0, 0, H, H I. Mokhun ct 
al use another definition of the Poynting vector singular- 
ity, according to which transversal component vanishes. 
This definition can be preferable, because Poynting vec- 
tor does not achieve exact zero at any point of the field 
.5] . It is evident that Mokhun's definition of the singu- 
larity includes our definition. In other words, we are to 
study only for a part of singular points which appear in 
terms of transverse component singularity. 

The most well-known type of singularity is the phase 
one. It arises, when the points of zero Poynting vec- 
tor are at the same time the points, where electric or 
magnetic field (or both) is zero, so that the field phase 
cannot be defined It seems to be the typical sin- 
gularity for scalar beams and finds a use for super- 
resolution 0, [^, [Tel . or construction and propagation 
of optical vorticesjl, Another well-known type of 
sing ularity is the polarization singularity (see theoretical 
01 El Q [3 and experimental [l6| results), when the 
state of polarization (in general, elliptical) is not fully de- 
termined. These singularities include the points of linear 
polarization, or L-points (rotation direction of polariza- 
tion is not defined) , and the points of circular polariza- 
tion, or C-points (polarization azimuth is not defined). 
At L- and C-points the Poynting vector is not zero in 
general. In some special situations, the energy flux den- 
sity may vanish owing to the certain directions of the 
vectors of electric and magnetic fields. It can happen at 
the points of linear or circular polarization. However, it is 



not the case in general, because vanishing of the Poynting 
vector due to polarization (polarization induced Poynting 
singularity) can be at the points of elliptical polarization, 
too. Finally, it should be noted that the fundamental 
connection between polarization and Poynting singulari- 
ties was noted in Refs. d @. Both singularities do not 
necessarily coincide, but go together. Such a relation can 
be applied for estimation of the magnitude of the angular 
momentum in some spatial region using field polarization 
characteristics. 

In Ref. [l| , the connection between Poynting singulari- 
ties and singular points of the dynamic systems has been 
developed. The main definitions and theoretical meth- 
ods related to arbitrary electromagnetic fields have been 
introduced therein. In the present paper we study the 
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singularities of vector Bessel beams. Bessel beams [17 1 
have been studied for several decades due to the prop- 
erties of non-diffraction, self-reconstruction, and angular 
momentum transfer (see, e.g., [H [H, [2(| and reviews 
[2lL |22| and references therein). In the current paper 
we consider the exact solutions of the Maxwell equations 
called vector Bessel beams. The fields of vector Bessel 
beams are characterized by some polarization distribu- 
tion, while the intensity of the beam has the form de- 
picted in Fig. [TJ A conventional Bessel beam follows 
from the vectorial one in the paraxial approximation. 
Some peculiarities connected with the vector nature of 
electromagnetic Bessel beams were discussed in the pa- 
pers [H [24|, HI [26[ . The number of applications which 
use both propagating and evanescent Bessel beams (such 
as super-resolution lenses mentioned above) grows from 
year to year. That is why it seems to be very important 
to study the regularities of spatial structure of complex 
electromagnetic fields in details. Singular peculiarities of 
electromagnetic beams (especially vortices) can be used 
as well, for example, in optical tweezers technique [27j |. 
detection of astronomical objects [28| . quantum cryptog- 
raphy [29], contrast enhancement in microscopy [3Q ]. etc. 

The special discussion is necessary on how to create a 
vector Bessel beam. The field with such intensity and po- 
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FIG. 1: Intensity distribution (longitudinal Poynting vector component in Eq. (|10[l ) of the TE-polarized vector Bessel beam 
with (a) m = 1 and (b) m = 2. Parameters: q/k = 0.5 (k is the wavenumber in vacuum), e — 1, and [i = 1 (dielectric 
permittivity and magnetic permeability of the medium). 



larization distributions arises in the core of the ordinary 
circular fibers. If the core radius is great, the intensity 
distribution contains a lot of minima and maxima. The 
field getting out the fiber can be considered as the vec- 
tor Bessel beam. The transverse wavenumber q of such 
a beam is defined by the dispersion equation. 

The paper consists of introduction, two sections and 
conclusion. Section 2 is devoted to propagating vector 
Bessel beams: we give the basic information about their 
origin from Maxwell's equations and consider the condi- 
tions of singularity generation. Section 3 deals with the 
evanescent beams which may arise under the conditions 
of total internal reflection and attenuate when propagat- 
ing. Both sections 2 and 3 include the discussion of elec- 
tromagnetic field properties at singular points (in partic- 
ular, the state of polarization) and the identification of 
the singularity type. 



II. SINGULARITIES OF PROPAGATING 
BESSEL BEAMS 

The m-th order Bessel beam solution of Maxwell's 
equations can be written in cylindrical coordinates (r, 

ip, z) as 



EM) 



H(r,iW lH(r, 



E(r,( 



exp(i/3z + im(f) — itut), (1) 



where longitudinal field components (along beam's axis) 
satisfy the Bessel equation 

£(£Ki(£)+<< 2 -^(£)=^> 

Here to is the circular frequency of electromagnetic wave, 
(3 is the longitudinal wavenumber, q = y/ k 2 efi — (3 2 is the 



transverse (radial) wavenumber, m is an integer number. 
Assuming finite- valued solution of equation ([2]) , it can be 
written as follows: 



H z = ciJ m (qr), E z = c 2 J m (qr), 



(3) 



where c\ and C2 are arbitrary complex numbers. Using 
these expressions for the longitudinal components, one 
can derive the other components of the fields [23, |2J] : 

E r = 5— cxJ m (qr) +i— c 2 {J m -x{qr) - J m+ i(qr)), 

q z r 2q 

E 4> = -' l Tr c i( J m-i{qr) - J m +i(qr)) - ^-c 2 J m (qr), 
2q q z r 

H r = i^-ci(J m _ x (gr) - J m+ i(qr)) + — ^— c 2 J m {qr), 
2q q z r 

H<t> = o~Ci J m (qr) + i— c 2 {J m -i(qr) - J m+1 (gr)J(.4) 

q z r 2q 

where k — lu/c is the wavenumber in vacuum, c is the 
speed of light. The meaning of the parameters c\ and c 2 
becomes clear, if we turn one of them to zero. If ci = 0, 
then H z = and the general expressions ^ and (0| are 
reduced to the fields of TM-polarized beam. Otherwise 
(c2 = 0) TE-polarized Bessel beam is obtained. It is very 
important that the amplitudes c± and c 2 are complex 
numbers which describe all possible vector Bessel beams 
of the order m. In contrast to commonly used scalar 
beams, one of the most important properties of these 
vector (electromagnetic) beams is polarization. 

The information to be discussed below is connected 
with the energy flux density (Poynting vector) of the 
beam considered. By substituting field components ([3]) 
and (g]) into S = (c/87r)Re(E x H*), the quantity S can 
be derived as 
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\m(ciC 2 )J' m J m ) e 



(5) 



where the derivative is calculated with respect to the en- 
tire function argument, i.e. J' m {qr) — dJ m /d(qr). This 
expression describes the dependence of Poynting vector 
on radial coordinate r for vector Bessel beam with any 
particular value of m which is a discrete parameter. This 
means that if the condition to = is taken (zeroth order 
beam), all terms with to in numerator vanish and there 
is no any ambiguities (and discontinuities) at the point 
r = 0. 

Two situations, m = and m ^ 0, are different from 
the mathematical point of view. Therefore, we will con- 
sider them separately. In the case to = 0, the expression 
for Poynting vector ([5]) is reduced to 

S r (r) = 0, 

S<t,(r) = -^-lm(c 1 c*)J (qr)J 1 (qr), (6) 

S,(r) = ^^\c 1 \ 2 +s\c 2 \ 2 )J 2 (qr). 

It is seen that, though the beam as a whole is directed 
along z-axis, the Poynting vector is not. This is due to 
vector nature of the beam considered which is a general 
solution of Maxwell's equations. Taking only TE- or TM- 
component (ci or c 2 is null) and proceeding to paraxial 
limit, one can obtain the usual (scalar) result with the 
only S z component. Note, that in Eqs. (JSJ) the first 
expression stands for the property of diffractionless of 
the Bessel beam, while the second one corresponds to 
angular momentum transfer by the beam. Finally, it is 
worth to stress that the straight physical sense of energy 
has not the Poynting vector itself, but the integral of it 



over the whole cross section of the beam. This integral is 
always directed along propagation of the beam as a whole 
(say, z-axis). But the Poynting vector of a complex beam 
(which consists of enormous number of plane waves) at a 
certain point can be directed in various ways (including 
opposite direction, —z) [23L [3ll| . 

The position of an optical Poynting singularity ro can 
be found from equation S(i"o) = or, in the case of prop- 
agating vector Bessel beam, 

S4r ) = 0, S z (r ) = 0. (7) 

The couple of equations ([7]) is satisfied, when 

/(ro) = Mgr ) = 0. (8) 

The positions of singularities ro are defined by the sin- 
gle beam parameter q. For paraxial beams (transverse 
wavenumber q is small compared with the wavenumber), 
the radial coordinates ro are much greater than that for 
non-paraxial vector beams. It is important that the sin- 
gularities exist for any c\ and c 2 . This feature is em- 
phasized in contrast to what will be discussed below for 
to ^ 0. 

If to 7^ 0, then using the formulae 

J'miqr) = \ (Jm-i(qr) - J m +i{qr)) , 

Jm ~Z (Jrn—1 "b Jrn+l ) 7 (9) 

2m 

the Poynting vector components (O can be represented 
in the form 



S r (r) = 0, 



cr 



32to,7t 



[Hli\ Cl \ 2 + e|c 2 | 2 )(J„ l _ 1 + J m+l f - 2/3Im(c lC ;)(J 2 l _ 1 - J 2 1+1 )] , 



^\ Cl \ 2 + e\c 2 \ 2 )(Jl_ 1 + Ji +1 )- 



1 + Pen 
2q~ 2 



Im(ciC2)(J, 



rn — 1 



-J, 



m+lj 



It leads to the two systems of equations (with upper and lower signs): 

/ ( r o) = (~/3± k^/eJI)J m -i{ra) + (/? ± k^/eJI)J m+ i(r ) = 0, 
H |ci| 2 + e|c 2 | 2 T2 A /£A*Im(ciC2) = 0. 

I 



(10) 



(11) 
(12) 



The solutions of equation pip are ro at different /3 > and m. Positions of singular points are determined by the 
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only coordinate ro in three-dimensional space. Therefore, 
the singularities are located on the cylindrical surfaces of 
certain radius. Reduced to one dimension, the surface 
can be considered as the point in the radial direction. 
The quantity of these surfaces is infinite due to infinite 
number of solutions of equation (|lip . The choice of sign 
in equation (jlip is connected with the choice of sign in 
relation (|12p which defines the coefficients c\ and C2 and, 
hence, the Bessel beam itself. Equation (fT2f can be rep- 
resented as quadratic form due to the common represen- 
tation of complex numbers ci = a\e 1Vl and c 2 = a^e^ 2 
(oi and a 2 are real and positive): 



\ia\ 



where Aip = ip\ 
Aip, we have 



ea 2 =F 2y / e/7etia2 sin Aip = 0, (13) 
Expressing the phase difference 



V2- 



sinA^ = ± 



SCLn 



(14) 



Taking into account that for oi, <22 > the inequality 
[ia\ + ea\ > 2^fejlaia2 holds, it should be stated that 
I sin Aip\ > 1. Owing to the limitation on value of sine 
equation (Ti"4"|) can be reduced to sin Aip = ±1, so that 



a 2 . 



Aip = ±'- 



2irp, 



(15) 



where p is integer. Thus, we conclude that vector 
Bessel beam can contain the singularities in beam's 
cross-section, only if it is constructed of TE- and TM- 
components oscillating out of phase and possessing am- 
plitudes matched by the wave impedance yj n/e. The 
simplest case of the relation (TT5]) (for propagation in vac- 
uum) gives C\ = 1, C2 = — i and c\ = 1, C2 = i for 
upper-sign and lower-sign equations (fl~5|) , respectively. 

As examples, we consider vector Bessel beams of the 
orders m = and m = 1. In figure Owe graphically solve 
equation ([5]) for m = and equation (jlip for m = 1. 
The zero crossings of the function /(r) (left-hand sides 
of equations ((8]) and (fTTj) ) in the upper figures corre- 
spond to the vertical arrows in the bottom figures. In 
figure (Hb) the most of Poynting singularities for the am- 
plitudes C2 = i and ci = — i appear together. The single 
exception is the first singular point for C2 = i. Spatial 
distance between a couple of Poynting singularities di- 
minishes with increasing the radial coordinate r, because 
the components of the energy flux densities are similar for 
great r. In figure[2fa), /(r)-curves for C2 = i and c 2 = — i 
(as well as for any amplitudes ci and C2) coincide. 

Since the radial component of the Poynting vector S r 
is absent (see equation (fl"U]) ), the direction of vector S 
in the plane {ip, z) can be specified by the angle ip{r) — 
axctan[S^(r)/S z (r)]. 

Figure |3fa) shows the change of the angle i[>(' r ) with 
radial coordinate. Its value oscillates between — tt and tt, 
the function tp being discontinuous. The function jumps 
from tt to — tt. For m = 0, the positions of jumps r\ 
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FIG. 2: Function /(r) (in the upper figures) and the compo- 
nents of Poynting vector (in the bottom figures) for (a) m = 
and (b) m = 1. Function / is defined as the left-hand side of 
equation (|8]) for m = or equation (jlip for m = 1. Selecting 
parameter c\ = 1 for m = 1, the another parameter C2 = =Fi 
follows from the choice of sign in equation Vertical ar- 

rows in the bottom figures mark the radii of singularities. 
Other parameters: (3 = 0.5fc, e = 1, fi = 1. 



coincide with that of singularities. In fact, at the singu- 
larities the expression for ip( r ) contains an ambiguity of 
the 0/0 type, which can be unwound as 



-0M = arctan^(r )/S^(r ), 



(16) 



where S'^ro) and S' z (ro) are the derivatives of the Poynt- 
ing vector components at the singular point ro. For 
in — we get 



47T 



/3Im(ciC2)Jo(<?r ) 



S' z (r ) = 0, 



so that ip{ r o) equals tt, i.e. singularity points coincide 
with points of jumps Tj. At the singularity only one of 
components of Poynting vector (namely, S^) changes the 
sign. Such a situation can occur, if S z is tangent to the 
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FIG. 3: Angle ij) vs. radial coordinate r for (a) m = and (b) 
m = 1. Vertical dashed lines mark the radii of singularities. 
Parameters: j3 = 0.5fc, e = 1, /x = 1, ci = 1, cz = — i. 



abscissa axis at singularity. This statement is confirmed 
by the figure Ula) for the Poynting vector components. 

The function ip(r) for m = 1 is shown in figure EJb). 
For some r, the jumps from —tt to tt occur. The posi- 
tions of jumps Tj are determined by vanishing S z (at the 
same time ^ 0). The radial coordinates of singular- 
ities Tq (see vertical dashed lines in figure [3Jb)) do not 
coincide with rj. It should be noted that at the points 
ro both components of Poynting vector, S z and S$, si- 
multaneously reverse sign. In figure [3|b) we observe that 
tp(ro) < 0, therefore, if 5^ > and S z < on the left of 
the singularity, then S$ < and S z > on the right of 
it. The distance between ro and rj gets smaller as radius 
increases, whereas the angle ip converges to — it (or 7r for 
another value of to). 

It should be noted that considered solutions for the 
vector Bessel beams are the centrally symmetric beams. 
They are very attractive due to the existence of the 
closed-form expressions for the fields. On the other hand, 
the studied beams are very specific, because of different 
imperfections during their realistic generation. So, if we 
introduce a small electric field e(r) as a fluctuation, the 
field can be written as E = E# + e and H = + h, 
where E b and H# are the electric and magnetic fields of 
the ideal vector Bessel beam, and h = l/(ifc/i)V x e. The 
fluctuation deforms the beam form, so that it becomes 
asymmetric. Then the Poynting vector takes the form 



S(r)«S fl (r) + s(r), 



4tt 



Re (Eg x h + e x Hg) 



It is evident that after adding s the singular point ro will 
be destroyed. Since the fluctuation field is an arbitrary 
one, it is quite impossible to meet the condition S(r) =0 
now (only for some specific fields e it can be done). How- 
ever, as the field e(r) is small, the singularities appear to 
be approximately situated at ro, however S(r) =0 is not 
satisfied exactly at any point. 



FIG. 4: Polarization parameter 7 versus radial coordinate 
for (a) m = 0, C2 = — i, (b) m = 1. The inset depicts the 
polarization ellipses at the points of solid curve marked by 
vertical arrows. Parameters: ci = 1, j3 = 0.5k, e = 1, (X = 1. 



Further we consider the electric and magnetic fields at 
the singular surfaces. None of them vanish as it would 
be in the case of phase singularity. In the case to = 
the azimuthal components vanish and both electric and 
magnetic fields are directed along e z , while c\ and C2 are 
arbitrary complex numbers, i.e. 

E(r, 0) = c 2 J Q {qr )e z , H(r, 0) = c 1 J (qr )e z . (18) 

Such longitudinal fields arc linearly polarized and can be 
treated as polarization L-singularities. 

For the rest of tos, one can derive the fields at singular 
points from equations ([3]) and 



E(r ,4>) = -c 2 J m+1 {qr )e 4 



H(r 0> < 



q 

+ c 2 J m (qr )e z , 
± ky/sjl 



-Ci- 



q 



Jm+iiqro^tj, 



+ C!J m (qr )e z 



(19) 



We may note that the magnetic and electric field 
strengths are connected as 



H 



£1 

C-2 



E. 



(20) 



By substituting the coefficients c\ and ci from equation 
(fT5)) . we represent the relation (pZU)) in the form H = 

Thus, the singularities considered are due to the spe- 
cial link between the electric and magnetic fields. We call 
them the polarization induced Poynting singularities. In 
general, they are distinct from the polarization singu- 
larities, which arise at the points of linear (L-points) or 
circular (C-points) polarizations [l3l [l4T|. 



The state of polarization can be studied by means of 
parameter 7 [32|: 



7 



\W 
IEI- 



(21) 



The introduced parameter characterizes the polarization 
of an electromagnetic beam using the complex electric 
field vector. The same analysis can be performed for the 
polarization of the magnetic field. All the types of polar- 
ization can be described with 7. So, circular polarization 
corresponds to 7 — 0, elliptical one — to < 7 < 1, 
linear polarization — to 7 = 1 . In figure [4] the singu- 
larities are located at the maxima of curves 7(7*) . That 
is why the electromagnetic fields at the singular surfaces 
are linearly polarized (7 = 1). The prove of 7 = 1 can be 
fulfilled analytically: the value 7 for the fields ((TSJ) equals 
I c 2 1 / 1 c 2 1 2 = 1. The inset of the figure SJb) clearly demon- 
strates the evolution of the forms of polarization ellipses. 
Starting with circular polarization at beam's axis (r = 0) , 
it changes step-by-step through elliptical to linear one at 
singularity tq. It should be noted that circular polar- 
ization appears only at beam's axis. Hence the axis is 
the place of C-points. Moreover, at r = the fields can 
be equal to zero (e.g., for to — 2), i.e. the phase singu- 
larity exists at this point. So, the beam axis simultane- 
ously corresponds to phase, polarization, and Poynting 
singularities. In this situation, the Poynting singularity 
is caused by the field vanishing, but not by the special 
polarization. However, the other Poynting singularities 
of the vector Bessel beam are the polarization induced 
ones. Sometimes, they can be identified as L-points ac- 
cording to the classification of polarization singularities. 
In our example, the Poynting singularities correspond to 
the L-points. In general, such is not the case, because 
the definitions of both singularities are not connected di- 
rectly. For m = the beam axis ro = is also the 



singular line (see figure |4]a,) , the electric and magnetic 
fields being non-zero ones and the Poynting singularity 
being caused by the linear field polarization. The link 
between polarization and Poynting singularities for arbi- 
trary electromagnetic fields is discussed in Refs. pi El- 



III. POYNTING SINGULARITIES FOR 
EVANESCENT BEAMS 

When a vector Bessel beam falls onto the interface be- 
tween two media, the total internal reflection can oc- 
cur. The Bessel beam that penetrates the interface ex- 
ponentially decays, when moving from the boundary. 
Such a wave is called the evanescent (inhomogeneous) 
one. Thus, the evanescent Bessel beam corresponds to 
the situation when q > k^JsJl. Applying replacement 

/3 -> i/3' = ivV - k 2 E[i to the fields © and @ we find 
the Poynting vector components [26j : 



S,p 

s z 



-^MdcS)^(«r)e-^-, 
C ^Mc^ + s\c 2 f)Jl( q r)e-"^ 



8nq 2 r 
cm(f3 12 + k 2 efi) 



(22) 

lm( Cl c;)J m (qr)J^(qr)e-^' z . 



In contrast to the propagating beam (see equation ©), 
the evanescent Bessel beam has the radial component S r , 
i.e. it is diffracted. As in previous section, one can mark 
out two cases, to = and to 7^ 0. 

For rn = 0, it is turned out that the Poynting vector S 
vanishes for each radial coordinate r. However, neither 
electric nor magnetic field is equal to zero. The fields are 
expressed as 



J 



E(r) = e _/3 z [ —c 2 .h(qr)e r + i — c t Ji(gr)e + c 2 J (qr)e z ) , 

1 q q 



H(r) = e 



, / Ql £ f c \ 

~ P 2 — ci Ji(qr)e r - i — c 2 Ji(qr)e^, + aJ {qr)e z 

\q q J 



(23) 



We have obtained the very exciting situation. There 
exist non-zero cylindrically symmetric solutions of the 
Maxwell equations, for which the energy is not trans- 
ferred (the energy flux is entirely zero) . Such a situation 
is similar to the case of the plane-wave standing waves. 
Being the Beltrami-fields, these standing waves do not 
transfer the energy, too [331 ]. 

Now the singularities are not located at certain cylin- 
drical surfaces. They occupy the whole three-dimensional 
space. So, the singularity can be not only the isolated 



point, or line, or surface, but the space itself. In general, 
the polarization of electromagnetic field (|23|) is elliptical. 
Exactly such a field polarization provides the complete 
vanishing of the energy flux. At the same time, polariza- 
tion is linearly or circularly polarized at some points. For 
example, L-points arise, when the condition J\ (qr) = 
holds true. It is obvious, that Poynting singularities in- 
clude the polarization singularities, which are situated at 
the cylindrical surfaces. 

For to 0, the Poynting vector components (12"2")) can 



7 



be rewritten using the properties of Bessel functions ([9]): 

S r = S^-I™(ciC*2)(Jm-l+Jm+l) 2 e^' Z , (24) 



16nm 
ckr 



327TTO 



(/x| Cl | 2 + e\c2\ 2 )(J m -i + J m+ ife- 2ti ' z , 



The singularity position ro is determined from the equa- 
tion 



f(r ) = J m -i(qr ) + J m +\{qr Q ) = 0. 



(25) 



The above equation holds for any complex coefficients C\ 
and c 2 . This means that any evanescent Bessel beam has 
singular points or, rather, cylindrical surfaces. In figure 
[SJa) the left-hand side of equation (f2"5| is shown. For 
m > I the Poynting singularity at beam's axis arises. 
It is caused by zero fields and coincides with the phase 
singularity as for the propagating electromagnetic Bessel 
beams. 

As to the electric and magnetic fields at singular 
points, they take the form 



E(r) = e im W' z J m+1 (qr ) ( ^c 2 e r + i f — Cl e, 

q q 



H(r) 



Armp—fi'z 



ft' 



ek 



J m+ i(qr ) —cie r - i — c 2 e,p . (26) 



q q 

It is important that the fields satisfy the relation 

H = aiE + a 2 E*, (27) 

where ot\ and a 2 are real and complex numbers, respec- 
tively. The numbers are equal to 



y\d\ 2 - e\c 2 \ 2 
2^Re(cic^) 



a 2 



lief + ec| 
2(iRe(cici,) 



Relation between the electric and magnetic fields (|27l) is 
necessary to turn the Poynting vector to zero. Gener- 
ally, the field link (|2"7| corresponds to the elliptical po- 
larization. The state of polarization can be described by 
the parameter 7 (see its definition (|2"T|) ). For evanescent 
Bessel beams, we will express it in the closed form as 



7 = 



|/? ,2 c 2 -^ 2 fc 2 c 2 | 
M 2 fc 2 | Cl | 2 + /3' 2 |c 2 | 



= / 2 A , 2 fc 2 /3' 2 [|c 1 | 2 |c 2 | 2 +Re(c 2 cf)] 

" V (M 2 fc 2 |ci| 2 +/3' 2 |c 2 | 2 ) 2 ■ 1 ' 

Figure [5^b) demonstrates three possible polarizations at 
the singularity point ro indicated by the vertical arrows. 
The linear polarization (L-point) is achieved for 7=1, 
i.e. for 



| Cl | 2 | C2 | 2 +Re( C 2 c^)=0 



(29) 
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FIG. 5: (a) Function /(r) (left-hand side of equation ([25|0 for 
different values of beam order m. (b) Polarization parameter 
7 versus radial coordinate for different amplitudes ci and C2, 
and m = 1. Arrows stand for singularity radii ro. Parameters: 
0' = 0.5k, s = 1, fi = 1. 



Using the above introduced representation of the complex 
amplitudes c\ = a\e llfl and c 2 = a 2 e 1V2 one can easily 
obtain the phase difference 

Atp = if! - <p 2 = | + 7rp, p = 0,±l,±2.... (30) 

In figure OJb) this condition is fulfilled for the amplitudes 
ci = 1 and c 2 = — i. 

Circularly polarized fields at singularities (C-points) 
are described by 7 = 0. From equation (|2"5)) it follows 
the necessary condition to obtain C-point at the Poynting 
singularity: 



2 c\ - n 2 k 2 cj = 0. 



(31) 



By substituting the amplitudes, we get 

Jf (32) 
Since ai and a 2 are real, the condition exp (2iA(^) = I 



8 



should be valid. Hence, 

A(p = np, p = 0,±l,±2... (33) 

and the relation between coefficients c\ and ci should 
satisfy the expression 

Such C-singularities are shown in figure O^b) as that cor- 
responding to the curve with amplitudes c\ = 1 and 
c 2 = 2. 

If neither equation (|3"0|) nor (j3~4")) holds, the polarization 
of the field at the singular point is elliptical. For exam- 
ple, the elliptically polarized field with 7 = 0.6 appears 
for ci = C2 = 1 (see figure Efb)). We can be convinced 
again that the Poynting singularity for the vector Bessel 
beams can include the polarization singularities as partic- 
ular cases. Of course, the polarization singularities can 
be generated at the points different from the Poynting 
singularities, too. 



strate the occurrence of non-phase optical singularities, 
namely Poynting singularities, which are situated at 
cylindrical surfaces in three-dimensional space. To cre- 
ate the singularity of the Poynting vector, TE and TM 
components of propagating Bessel beams defined by the 
complex amplitudes c\ and ci should be out of phase and 
be matched with the wave impedance. Such a strong re- 
striction on the amplitudes arises only for propagating 
waves. Any evanescent Bessel beam is appeared to be 
singular. Moreover, the evanescent beam with m = is 
singular at any spatial point. This case describes a sort 
of cylindrical standing wave. The points of Poynting sin- 
gularities always include the phase singularities as partic- 
ular case (at Bessel beam's axis). However this is not the 
case for L- and C-points (polarization singularities). The 
manifold of the Poynting singularities and manifold of L- 
and C-points can intersect, however, it is only the coinci- 
dence. Since vector cylindrical solutions of the Maxwell 
equations describe the electromagnetic modes of circular 
fibers, the theory of Poynting singularities can be applied 
for the fields of guiding structures. 



IV. CONCLUSION 



Summing up, vector Bessel beams (solutions of 
Maxwell's equations in cylindrical coordinates) demon- 
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